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Density-density functionals and effective potentials in many-body electronic structure calculations 
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We demonstrate the existence of different density-density functionals designed to retain selected properties 
of the many-body ground state in a non-interacting solution starting from the standard density functional theory 
ground state. We focus on diffusion quantum Monte Carlo applications that require trial wave functions with 
optimal Fermion nodes. The theory is extensible and can be used to understand current practices in several 
electronic structure methods within a generalized density functional framework. The theory justifies and stimu- 
lates the search of optimal empirical density functionals and effective potentials for accurate calculations of the 
properties of real materials, but also cautions on the limits of their applicability. The concepts are tested and 
validated with a near-analytic model. 
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I. INTRODUCTION 

Density Functional Theory (DFT)^ is based on the 
Hohenberg-Kohn proof of a functional correspondence be- 
tween the ground state energy and the ground state density 
E[p(r)]. In the formulation of Kohn and Sham (K-S))S the in- 
teracting electron gas is replaced by non-interacting electrons 
moving in an effective potential. In this construction, the non- 
interacting density p(r) is equal to the interacting one p(r) but 
no other property of the interacting ground state is in principle 
retained in the non-interacting wave function. Initially DFT 
was formulated to describe the total ground state energy of an 
interacting system and p(r)w^£ Although progress towards 
more accurate density functionals is ongoing, current approxi- 
mations such as the local density approximation (LDA)^ and 
more recent gradient-based extensions^ are already success- 
ful in predicting many electronic properties of real materials. 
This success has led to the use of DFT beyond its formal scope 
and unfortunately tempted some to believe that if we had the 
exact ground state density functional, we would only need to 
solve non-interacting problems even for properties not related 
to the ground state energy and density. While the virtues and 
limitations of the Kohn-Sham eigenvalues are discussed in 
textbooks^ 6 - the possible reasons for the success or failure of 
Kohn-Sham wave-functions in many-body problems are little 
understood and widely dispersed throughout the literature. 

It is often assumed, without a formal proof, that the Kohn- 
Sham non-interacting ground state wave-function forms a 
good description of the interacting ground state wave-function 
to be used as the foundation of theories that go beyond DFT 
such as GW-Bethe Salpeter Equation (GW-BSE), Quantum 
Monte Carlo (QMC), or even configuration interaction (CI). 
This leads to an apparent contradiction in the literature since 
density functionals that provide wave-functions that are a 
good starting points in one field (as judged by comparison 
with experiment) are found inadequate in others. Broadly 
summarizing: for structural properties gradient corrected den- 
sity functional* 7 - are nowadays preferred over LDA. 2,4 Struc- 
tural properties depend essentially of atomic forces which in 
turn are related to the density. However, for GW-BSE cal- 



culations of optical properties, an LDA-based ground state 
is preferred^. In this approach a good initial approximation 
for the Green function is required. In QMC calculations a 
(non-interacting) Hartree-Fock (HF) ground state might be 
preferred over LDA, but the subject is still under debate. In 
CI calculations, instead, it is empirically claimed that the con- 
vergence with natural orbitals^ is more rapid than HF orbitals. 

In Diffusion Quantum Monte Carlo (DMC) a trial wave 
function enforces the antisymmetry of the electronic many- 
body wave function^ and the nodal structure of the solution. 
The accuracy of the trial wave-function is critical and deter- 
mines the success or failure of the method to accurately pre- 
dict properties of real materials. The trial wave-function is 
usually a product of a Slater determinant $r (R) and a Jastrow 
factor e J (R). <j>r(R) is often constructed with single parti- 
cle Kohn-Sham orbitals or from other mean field approaches 
such as HF. The Jastrow, in turn, is a symmetric factor which 
does not change the nodes, but accelerates convergence and 
improves the algorithm's numerical stability. The DMC algo- 
rithm finds the lowest energy of the set of all wave-functions 
that share the nodes of ^y(R). The exact ground state en- 
ergy is obtained only if the exact nodes are provided. Since 
any change to an antisymmetric wave-function must result 
in a higher energy than the antisymmetric ground state, the 
energy obtained with arbitrary nodes is an upper bound to 
the exact ground state energy^ Only in small systems is it 
possible to improve the node a 11,12 ' 13 ' 14 ' 15 or even avoid the 
trial wave-function approach altogether 16,17 . Consequently, 
a general formalism that could alleviate the nodal error in 
large sy stern a 18,19 is highly desired. Quite recently it has been 
shown that within the single Slater determinant approach the 
computational cost of the DMC algorithm can have an al- 
most linear scaling with the number of electron a 18,19 ' 20 ' 21 . It is 
claimed, if not formally proved, that the nodes of the many- 
body wave-function are not too far from those of a wave func- 
tion obtained via a mean field approach. However, this might 
not continue to hold as electron-electron interactions become 
more important. To improve the accuracy of these approaches 
and increase the range of materials to which they can be ap- 
plied it is important to examine the advantages of different 
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mean-field wave-functions. 

In this paper, we demonstrate that density-density function- 
al can be obtained by finding the minimum of different cost 
functions relating the set of non-interacting w-representable 
ground state with an interacting many-body state. The mini- 
mum of these cost functions establishes a correspondence be- 
tween the non-interacting and the interacting wave-functions 
and their associated densities and potentials. The cost func- 
tion can be designed to retain selected properties of the many- 
body wave function in the non-interacting one. Crucially, for 
DMC applications the nodes can be optimized. Under cer- 
tain conditions density-density functionals exist that can lead 
to standard scalar-density functionals. As in the case of stan- 
dard DFT, this proof does not mean that we know the expres- 
sion of each functional or associated potential but it will cer- 
tainly stimulate the search of methods to find or approximate 
them. For DMC applications it is enough to prove that an op- 
timal mean field potential for nodes exists. In order to test 
the theory, we find the ground state wave function of a model 
interacting system. Then we obtain (i) the exact DFT effec- 
tive exchange correlation potential associated with the ground 
state density, (ii) the potential that maximizes the projection 
of <1>t(R) with the ground state. Finally, (iii) we optimize 
a potential to match the nodes and find that surprisingly, for 
this model, the non-interacting solution in the same potential 
as the interacting problem is a very good approximation for 
the nodes while the exact non-interacting Kohn-Sham ground 
state is particularly poor. 

This paper is organized as follows: In Section II we demon- 
strate the existence of an different density-density correspon- 
dences associated with cost functions. We prove the exis- 
tence of this functional correspondence for the case of optimal 
nodes required in DMC. In Section III we solve an interact- 
ing problem up to numerical precision and find its many-body 
ground state wave-function. Subsequently we optimize differ- 
ent cost functions to retain specific properties of the ground 
state. Finally in Section IV we discuss the relevance of our 
results for many-body electronic structure and give our con- 
clusions. 



H. GENERALIZED DENSITY-DENSITY FUNCTIONALS 

Given an interaction in a many-body system, the 
Hohenberg-Kohn theorem^ establishes a functional cor- 
respondence between densities p(r), external potentials 
V(?")[j0(r)] and ground state wave-functions $(R)[p(r)]; 
where [p(r)] denotes a functional dependence on the ground 
state density, and R denotes a point in the many-body 37V 
space. Since the density changes according to the strength 
and functional form of the interaction, this correspondence 
is different for different interactions. For a fixed interac- 
tion, the subset of densities p(r) corresponding to a ground 
state of an interacting system under an external potential V(r) 
are denoted as pure state w-representable4 A non-interacting 
pure state u-representable density is given instead by p(r) = 
Si/ \$v (r) 1 2 where <j) v (r) are the Kohn-Sham-like single 



particle orbitals, or eigenvectors, of the Hamiltonian: 



-\v 2 + V{v) 



4> y (r) = e v <\> v (r) . 



(1) 



where V (r) is an effective single particle potential. For sim- 
plicity we denote here a density to be w-representable if it is 
both pure state non-interacting and pure state v-representable. 
In the following we also imply pure state when we write only 
v-representable. 

Each point in the sets of v-representable densities is asso- 
ciated with two different points in the wave-functions Hilbert 
space. In figureQ]we schematize the subset of u-representable 
densities and the functional correspondence with the sub- 
sets of the interacting and non-interacting ground state wave- 
functions. Note that, in principle, the two subsets of ground 
state wave-functions do not necessarily overlap. In the non in- 
teracting case the wave-function is given by a Slater determi- 
nant of Kohn-Sham-like orbitals but for interacting problems 
this simplification is not longer possible. 

The Kohn-Sham scheme for density functional theory es- 
tablishes a correspondence between interacting and non- 
interacting wave-functions represented as line (1) in Fig Q] 
This Kohn-Sham correspondence between wave-functions is 
implicit in the Khon-Sham construction for the external effec- 
tive potential £ Figure Q] emphasizes that the wave-functions 
joined by line (1), while different, give the same electronic 
density. In more technical terms they both belong to the 
same Percus-Levy partition of the Hilbert spac o 22 ' 23 but they 
are the minimum energy wave-function for different interac- 
tions. The exchange-correlation potential is by construction 
the difference that one has to add to the external potential in a 
non-interacting problem so that its ground state density is the 
same as the interacting one. If the energy-density functional 
E[p(r)] is known, the effective non-interacting potential can 
be obtained following the standard Kohn-Sham approach. If 
the ground state density p(r) is known, the same correspon- 
dence between interacting and non-interacting densities can 
be achieved by minimizing the following function 



P(r)Y 



(2) 



within the subset of non-interacting u-representable densities. 
Formally, this could be done by exploring all values V (r) in 
Eq©. 

In practice, if the density of the interacting ground state is 
known, the potential Vk (r) that minimizes Eq. (01 can be 
obtained numerically with a procedure similar in spirit to the 
optimized effective potential (OEP) method. The change in 
the density required to minimize Eq. (ffj) is 



\ = -[p{v)-p{r)\. 



(3) 



Within linear response, the change in the potential required to 
produce A p is 

AV Kp (r) = / dr' W) - p(v')] S -^l (4) 



6p(r) 

Adding recursively AVk p (r) we can find the potential Vk (r) 
associated with K p = (see an example below). 
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FIG. 1: (Color online) a) Representation of the sets of pure state 
w-representable interacting densities, b) Sets of interacting and non- 
interacting ground state wave-functions. The Kohn-Sham formula- 
tion of DFT relates au - representable density with a pair interact- 
ing and a non-interacting wave-function. The same functional corre- 
spondence can be obtained minimizing Eq. (O (line 1 in the figure). 
Different cost functions relate an interacting ij-representable density 
with a different non-interacting-^ representable density (see lines 2 
and 3). 



A. Other density-density correspondences 

It is often desirable o preserve properties in addition to 
the density of the many -body ground state ^(R) in a non- 
interacting wave-function <£>t (R) to be used as a starting point 
for theories that go beyond DFT. This task involves exploring 
all the non-interacting u-representable set in order to find a 
wave-function that best describes a given property. This is 
a typical optimization problem. One of the most common 
strategies in optimization is the design of a cost function. One 
example is Eq (|2), a measure of the difference in two densities. 
Another example of a cost function is 



A" 



Det 



(5) 



which involves a projection of the interacting ground state ^ 
in the set of non-interacting w-representable wave-functions 
{<£>t}. The minimum of Eq. © is the non-interacting ground 
state Slater determinant with maximum projection in the in- 
teracting ground state. We have claimed above that the inter- 
acting wave-functions might be in general very different from 
a single non-interacting Slater determinant. Accordingly, we 
expect Kn et > — 1. 

We expect to find a different minimum in the non- 
interacting ground state set, if we change the functional form 
of the cost function from Eq. (O to Eq. (© for the following 
reasons: 

1) We can visualize the cost function as a scalar potential 
defined in the full Hilbert-space. Although different cost func- 
tions can share the same minimum in the complete Hilbert 
space, in the restricted subset of non-interacting ground state 
wave-functions, different cost functions can have a different 
minimum: the optimal point found depends on the functional 
form of the cost function. Accordingly, while all the cost 
functions we propose here would be minimized if we could 



reach the interacting many body state "J ( where, of course, 
every property is retained exactly), because our search is con- 
strained to non-interacting u-representable subset the mini- 
mum we would find will depend on the properties we wish 
to retain. 

2) The Hohenberg Kohn theorem, when applied to the non- 
interacting u-representable case implies that, in the absence 
of degeneracy, there is at most one wave-function that has the 
same density as the interacting case. Therefore, once the min- 
imum of an arbitrary cost function is found, its associated 
non-interacting density can no longer be equal to the interact- 
ing density unless the property enforced by the cost function 
can be related back to density. Enforcing the non-interacting 
density to remain equal to the interacting ground state den- 
sity prevents all other properties of the non-interacting wave- 
function from being further improved. If we intend to opti- 
mize other properties, we have to relax the density constraint 
finding a different wave-function associated with a different 
density. 

The minimization of different cost-functions, relating 
the interacting ground state \I/(R) with the non-interacting 
u-representable set, provide in-principle different corre- 
spondences between interacting and non-interacting wave- 
functions represented as different lines in figure Q] Each 
cost function K defines a correspondence different than 
the identity between pure state w-representable densities and 
non-interacting pure-state v— representable densities. As a 
consequence, the idealized optimization processes outlined 
here defines an operator Uk that turns each p(r) into a 
non-interacting density corresponding to the wave-functions 
<I>t(R) which is the minimum of a cost function K. 



Pk(v) = U K [p(r)] . 



(6) 



Note that if the minimum of a given cost function K is a 
single $t (R) for every w-representable density, then Uk de- 
fines a density-density functional. When more than one non- 
interacting u-representable wave-function give the same opti- 
mal value for K, the degeneracy can be broken by additional 
requirements in the cost function [e.g. also minimizing Eq. 
(O, the difference between the current and pure state densities 
]. Since we only need one optimal wave-function, any from a 
degenerate minimum can be chosen to construct the density- 
density functional Uk- When minimization of a cost function 
defines a one to one correspondence with an inverse, a more 
usual energy-density functional of the form E{U^ [pk (r)]} 
can be constructed. Only a restricted class of cost functions 
lead to density transformations with an inverse. Minimization 
of the cost functions among all pure-state-non-interacting v- 
representable densities defines the optimal effective potential 
which is a function of this density. 

Given a cost function, K, finding an approximation for the 
density transformation operator Uk could certainly be as de- 
manding as finding an approximation for the energy-density 
functional E [p{r)] required by DFT based methods. This task 
is beyond the goal of this paper. However, we will show that 
we can expect the operator associated to the best nodes for 
DMC (Udmc) to be non-local and very different from the 
identity. Accordingly we can expect non-interacting wave- 



4 



functions with good nodes to be a poor source of densities. 
Moreover, for the example considered below, we find, that the 
direction we might have to explore to optimize the potential 
might be surprisingly different than the attempts considered 
so fa r 21 ' 24 . 



B. The Diffusion Monte Carlo case 

We next show that optimization of the nodes for DMC 
among the set of v— representable wave-functions leads to a 
correspondence between pure state v-representable densities 
and pure state non-interacting v-representable densities of the 
class described above. These in turn demonstrate the existence 
of an optimal effective non-interacting nodal potential. 

Since, the ground state density p(r) determines the ground 
state wave-function ^(R)[/9(r)]r. p(r) defines also the points 
R of the nodal surface Sb(R)[p(r)] where vp(R)[/j(r)] = 
0. We can also classify the nodal surfaces in pure state v- 
representable and pure-state-non-interacting u-representable. 

The DMC algorithm in the fixed node approximation finds 
the lowest energy of the set of all wave-functions that share 
the nodes or the trial wave-function. For Slater determinant 
Jastrow wave-functions, the nodes of the trial wave-function 
are by construction those of $^(R); that is they are pure-state 
non-interacting v-representable. The DMC energy, Eqmc is 
also a function of the external potential which in turn is a 
function of the interacting ground state density V(r)[p(r)]. 
Thus minimization of Ecmc[^t(R), P( r )] m tne set of 
non-interacting v-representable wave-functions $y(R) deter- 
mines one $t(R) with the best nodes. Every optimal $x(R) 
defines an optimal auxiliary density p~DMc( r )- As a conse- 
quence optimizing the nodes of the trial wave-function by per- 
turbing the nodes of pure state non-interacting wave-functions 
implies finding another correspondence between interacting 
and non-interacting densities (another line in figure [TJ. The 
best cost function for optimal nodes is ultimately the DMC 
energy itself. 

Since we restrict the search to pure-state non-interacting v- 
representable nodes, the minimum energy Edmc [p( r ] w iU t> e 
larger than the true ground state energy E[p(r], because of 
the upper bound theorem, unless So (R) is non-interacting v- 
representable. 

Note that for an arbitrary interaction Sq (R) is not expected 
to be, in general, pure-state-not-interacting v-representable. 
However, if Sq (R) were non-interacting v-representable, the 
best Slater Determinant $t(R) for DMC could be formally 
found by finding the minimum of the cost function 

K So = f dS|$ r (R)| 2 . (7) 

J S 

where J s denotes a surface integral over the interacting nodal 
surface. 



III. COST FUNCTION MINIMIZATION 

To demonstrate the theoretical concepts above we solve a 
simple non-trivial interacting model as a function of the in- 
teracting potential strength and shape. We then optimize the 
wave-functions to minimize the cost functions in Eqs. (0, 
(0 and (0 so as to find the exact DFT wave-function, the 
wave-functions that maximize the projection on the interact- 
ing ground state and minimize the projection on the nodes. 
Subsequently, we estimate the volume of the Hilbert space 
enclosed between the nodes of the interacting wave-functions 
and the optimized non-interacting ones. 

A. A model interacting ground state 

For illustrative purposes we choose the interacting prob- 
lem to be as simple as possible and yet not trivial. We 
solve the ground state of two spin-less electrons mov- 
ing in a two dimensional square of side length 1 with 
a repulsive interaction potential of the form F(r, r') = 
87 cos [aTr(x — x')] cos [an(y — y')].— While this potential is 
different than the Coulomb interaction, it shares some of its 
properties. For positive 7 and \a\ < 1 the interaction is repul- 
sive with a repulsion that increases monotonically when for 
shorter distances. The amplitude of the repulsion as compared 
to the kinetic energy can be changed by adjusting 7. Since the 
Coulomb interaction is self similar, changing 7 mimics what 
happens in a real system when we change the size of the sys- 
tem. The shape of the potential within the confined region can 
be altered by changing a. In the limit of a — > the interaction 
potential is separable which allows several limits to be tested 
(such as the nodes). The functional form facilitates an analyt- 
ical treatment of the problem by removing the singularity of 
the Coulomb interaction at short distances. 

We expanded the many-body wave-function in a full CI 
on non-interacting Slater determinants with the same sym- 
metry as the ground state. The ground state is degenerate 
because there are only two electrons. We chose one of the 
ground state wave-functions according to the D% subgroup 
of the D4 symmetry of the Hamiltonian. With this choice, 
p(r) has D2 symmetry (x is not equivalent to y). The basis 
of Slater determinants was constructed with functions of the 
form 16sin(n7rx)sin(m7r?/) with n and m < 8. Since parity is 
preserved by the interaction and Slater determinants of iden- 
tical functions are zero, the size of the basis set is reduced to 
only 300 in our calculations. 

Most of the calculations reported here were done analyt- 
ically with the help of the Mathematica package, including 
all of the electron-electron interaction integrals^ The only 
source of errors are numerical truncation and the size of the 
basis, which was tested for convergence. 

In Figure|2]we show the quadrants of densities correspond- 
ing to wave-functions that are even for reflections in the y 
direction and odd in the x direction for 7 = 2 and a = 1. 
Figure|2ja) shows the upper left quadrant of the density of the 
interacting ground state of two spin-less electrons obtained 
with full CI. Figure|2jb) shows the upper right quadrant of the 
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FIG. 2: Ground state densities, in particles per unit area, for two in- 
teracting spin-less electrons in a square box. The complete density 
can be obtained by reflections (see also Fig[3j. (a) Full CI ground 
state interacting density, (b) Exact DFT solution obtained minimiz- 
ing Eq. l(2}. (c) Slater determinant with maximum projection with CI 
the ground state [see Eq. J5j]. (d) Slater determinant with minimum 
amplitude on the nodes of the CI ground state [see Eq. 10]. 



non-interacting density corresponding to the effective poten- 
tial obtained by adding recursively AVk p (r) [Eq. ©], which 
is exactly the reflection of p(r') Fig. 0a) up to numerical pre- 
cision [K p = in Eq. (0]. Because of the Hohenberg-Kohn 
theorem^ Vpc„ and the Kohn-Sham potential Vks (p) can only 
differ by a constant and thus the wave-functions coming from 
this potential are the exact DFT wave-functions for our inter- 
action. The properties of the wave-functions will be discussed 
later in the text. The densities in Figs. 0d) and0d) corre- 
spond to the minimum of the cost functions given in Eqs 
and obtained as described below. 



B. Effective potential optimization 

We now consider more difficult cost functions than den- 
sity differences, Eq. 01. For non-interacting u-representable 
densities there are also functional correspondences between 
ground state wave-functions, potentials and densities. This 
concept has been exploited in the optimized effective poten- 
tial (OEP) for exact exchange ^ 27,28,29 ' 30 The exchange poten- 



tial can be calculated in OEP 



27.28.29.30 



as: 



V x (r) 



E 



5E X 
Sp(r) 



dr'dr" 



K (r") 



K (r") SVks (r') 



c.c. 



(8) 
SV KS (r') 



Sp(r) 



In Eq. (0, the functional derivative SE x jS(\) v (r) is evaluated 
directly from the explicit expression for the exchange energy 
E x in terms of <)>„ (r). Next 8<t> v (r) /SVks ( r ') is evaluated 
using first-order perturbation theory from Eq. (0. Finally 



SVks ( r ')/Sp (r) is the inverse of the linear susceptibility op- 
erator. If there are fixed boundary conditions such as the num- 
ber of particles, the susceptibility operator is singular^ 3 ^. Ex- 
cluding these null spaces it can be inverted numerically i2£ We 
use earlier this susceptibility operator in Eq. (0. Equation 
(0 is by construction the gradient of the exchange energy in 
the set of pure-state-non-interacting u-representable densities. 
While we are not going to attempt an exact exchange approach 
in this paper, the ability to calculate gradients allow us to min- 
imize cost functions as long as the cost function K can be ex- 
pressed in terms of non-interacting ground state wave wave- 
functions or eigenvalues. The potential that minimizes K can 
be obtained by recursively applying the formula 



8V K {v) 



OCC 



dr'-; 



SK 



k (r') 



(r') SVks (r) 



+ c.c. 



(9) 



Equation (0 gives the direction we need to change the po- 
tential to minimize the cost function. The magnitude of 
the change is controlled by e, which can be adjusted as one 
reaches the minimum. 

Replacing K by Kjj e t m Eq. (0 and using Eq (0 and first 
order perturbation theory we find 



SV Koct (r) = e(*|$ T )^^(*|ctc,|$ T ; 



</>n(r)</w(r) 



+ C.C. 



(10) 



In equation ([TUt J2n ( ) means sum over occupied (un- 
occupied) states, while c^ n and c v are creation and destruc- 
tion operators on the non-interacting ground state |$t)- One 
can understand also the state c^c^ |$t) as the many body 
wave-function resulting from replacing the occupied 

state <f) y by the <j> n . This is equivalent to creating an electron 
hole pair excitation in a non-interacting ground state. In Eq. 
( fTOb a term in the potential is added every time an electron 
hole pair excitation has no zero projection to the interacting 
ground state. Since the basis of products of wave-functions 
4> n {r)(f) u {r) is over-complete, there are linear combinations 
with non-zero coefficients that add up to zero. A minimum 
is found when the gradient of the cost function with respect 
to variations of the effective potential is zero. If the absolute 
minimum is found, the wave function can only be improved 
further by a multi-determinant expansion, that is, outside the 
set of pure-state non-interacting densities. Since we choose a 
basis expansion for the single particle orbitals to be sine func- 
tions, the products cf> n (r)(f> nu (r) are linear combinations of 
sine products. These sine products can be transformed analyt- 
ically to cosines. The change in the potential is thus written in 
a cosine basis which is complete. All coefficients must van- 
ish in the cosine basis when a minimum is found. This allows 
us to verify that the gradient in the potential can be minimized 
up to numerical precision. The integrated effective potential is 
thus naturally expressed as a linear combinations of cosines, 
which allows the analytical calculation of the coefficients of 
the effective potential matrix in a basis of sines, where the 
kinetic energy is diagonal. The Slater determinant |$y > is 
written in the same basis as the interacting ground state |$ >. 
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The projections involved in Eq. ([TUt are then reduced to a 
scalar product of the vectors of coefficients. 

In figure W[c) we show the ground state density associated 
to minimization of Eq. (0 for the same parameters as the in- 
teracting ground state density in Fig 0a). We see that while 
optimizing the cost function (0 allows matching the interact- 
ing density exactly, optimizing the wave-function projection 
requires a significant change in the resulting density. 

Similarly, replacing K by Ks a in Eq. (0 and using Eq (0 
we get 



j^jZ I dS$^(R)$ T (R) 



+ c.c. 



0n(r)0rm(r) 

(ID 



Unlike Eqs|4]and[T0] a complication appears when evaluating 
the integral over the nodal surface J s dS. This integral in- 
volves finding the points where the many-body wave-function 
is zero. The problem is simplified because the derivatives 
of the many body wave-functions can be obtained analyti- 
cally. Consequently, starting from an arbitrary point R we 
can find a zero recursively with the Newton-Raphson method, 
R Il+ i = R„ + *(R)W(R)/ |V*(R)| 2 . Next we make 
a random displacement AR in the hyper-plane perpendicu- 
lar to V4'(R) within a circle of radius 0.05 and find a node 
again. We repeat this process 50 times and select an element 
of {R}s- With this parameters, the random position R„ can 
travel across the full size of the system so that the distribution 
is homogeneous. By repeating this process N = 500 times, 
excluding points at the boundaries which are zero by construc- 
tion, we generate an homogeneous distribution of points at the 
nodal surface {R}s- We approximate the integral in Eq. 
as a sum on the values on the set {R}s- Note that while the 
total area of the surface would be in general involved as a fac- 
tor, the value of this area is not relevant since we are interested 
in finding a minimum of the cost function and the position of 
the minimum of any function is not altered by a positive mul- 
tiplicative constant. The sum over random points introduces a 
relative error of order 1 / V500 . Replacing the integral with a 
summation creates also many local minima in the landscape of 
Eq. ( fTTT i. Accordingly, we tested different initial conditions; 
the best results are obtained starting from Vk s = 0. 

The density resulting from minimization of Eq. (0 is plot- 
ted in Fig 0d). We see here again a significant change as 
compared with the fully interacting CI ground state [see Fig 
0a)] and the the exact DFT non-interacting solution Fig0b). 

Figure is a clear example that corroborates our claim 
in Section III Al that enforcing different properties on the 
non-interacting wave-function implies a density-density cor- 
respondence different than the identity between the interact- 
ing and non interacting systems. Similar results are observed 
as function of the strength and shape of the interaction [con- 
trolled by 7 and a] 

A comparison between Eqs. (TTOb and (fTTT i clearly 
shows that the relative values of the coefficient multiplying 
(f> n (r)(j) nu (r) depends fundamentally on the cost function. 
Therefore, even starting from the same effective potential 
and $t(R) the coefficient affecting each individual product 



4>n (r)</>nu ( r ) depends on the functional form of the cost func- 
tion. This change in the potential remains present when the 
potential is written in the complete cosine basis. Thus, the ef- 
fective potential must change in accord with the property of 
the interacting ground state that one aims to enforce in the 
non-interacting ground state with a cost function. 

Figure0shows the effective potentials used for the calcula- 
tions shown in Fig [2] We show in Fig 0a) a constant, since in 
the interacting problem solved with full CI no effective exter- 
nal potential was added. Figures 0b) [minimum of Eq. (0], 
13c) [minimum of Eq. (0] and [3jd) [minimum of Eq. (0] 
show a clear change in the effective potential depending on 
the cost functions. As argued earlier Fig. Ob) shows the exact 
Kohn-Sham DFT potential for this interaction which implies 
that a different density functional must be used to obtain non- 
interacting wave-functions preserving properties other than 
the density. 

Note that Eq. could be zero only for a pure-state-non- 
interacting v-representable nodal surface. However, if the 
nodes are not, replacing in Ks could result in a potential that 
simply prevents the non-interacting wave-function to reach re- 
gions of space where the nodes are more troublesome. The 
potential shown in Fig. 0d) presents a maximum in regions 
where instead Figs. 0b) and0c) develop a minimum. These 
are the regions where the electrons in the many-body wave- 
functions tend to localize because of correlation effects. The 
maximum in Fig [3d) suggest the possibility of non-u repre- 
sentability by a non-interacting wave-function in this model. 



C. Wave-function internal structure 

In order to quantitatively test the quality of the nodes of the 
wave-functions found by minimization of Eqs. (gj, (0 and 





^1 (b> 


I© i • ' 





FIG. 3: Optimized effective potentials corresponding to the densities 
in Fig. [2] The complete potentials can be obtained by reflection 
on the black lines. Gray level valued are given on the right. The 
optimal effective potentials are strongly dependent on the property 
we target to retain in the wave-function. 
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and to test the convergence of the nodes of the full CI calcula- 
tions, we take advantage of the homogeneous distribution of 
points {R} s at the nodal surface So (R) described earlier. 

For each point R in {R}s we can find the distance £i to 
the node of another wave-function $t(R) in the direction of 
W(R). Thus 



AV 



1 



(12) 



is an approximated measure of the fraction of the Hilbert 
space volume between the nodal surface of $^ (R) and ^ (R) 
and 



6p 



1 

3iV 



$>|3>r(Ri)f 



(13) 



measures the probability density inside AV. 

We can use Eqs. ( fT2l and ( fT3T > to test the convergence of 
the CI ground state nodes as a function of the size of the basis 
set. While the ground state energy requires 40 basis func- 
tions, the nodes are more difficult to converge requiring four 
times as many. The size of the basis required to converge the 
nodes was determined for 7 = 2 plotting 5p between the CI 
ground state with 300 wave-functions and the CI ground state 
obtained using a reduced basis. 

The quantities in Eqs. dl~2t and (fljjl can be used also to 
characterize the nodes of different wave-functions as com- 
pared with the exact node. In figure |4] we show the volume en- 
closed between the nodes AV of different optimized $t(R) 
and the interacting ground state ^(R) as a function of the 
strength of the interaction potential 7. Note that the Kohn- 
Sham DFT solution gives a significantly larger volume than 
other optimized wave-functions. The difference increases as 
the interaction strength increases. The wave-function that re- 
sults from minimizing Eq. (O which targets wave-function 
projection fares very well over the range explored. In turn, 
minimization of Ks [see Eq. (|7]i] results in nodes that 
are only sometimes marginally better. Surprisingly, the non- 
interacting solution, that is the non-interacting ground state 
in the absence of any effective potential, is remarkably good. 
Similar results are found by altering the shape of the potential 
with a. 

In figure [5] we plot the values 5p for different optimized 
wave-functions as a function of 7. We see again that the exact 
Kohn-Sham DFT solution is not the best. The quantity Sp is 
a measure on how much the error in the nodes would affect 
the probability density and thus it can be understood as the 
as a measure of the nodal error in the ground state energy. 
Again in this case the non-interacting ground state without 
any effective potential is the best approximation. 



IV. DISCUSSION 

Although the numerical investigation of the different 
density-density functionals described above required numer- 
ical representation of the many body ground state wave- 
function, the conclusions that we draw have general value. 



" 0.015 

4 



5 0.010 
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0.5 1 1.5 2 

Interaction strength 7 

FIG. 4: (Color online) Fraction of the Hilbert space A V between the 
full CI node and the nodes of different optimized wave-functions. 
Triangles correspond to the exact DFT wave-function [Eq. (2)], 
squares to maximum projection [Eq. (3)], rhombi to the minimum 
amplitude at the nodes [Eq. ([7J] and circles to the non-interacting 
ground state. The inset shows the method used to estimate AV" 




0.5 1 1 .5 
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FIG. 5: (Color online) Probability density inside the volume between 
the nodes of the full CI wave-function and optimized wave-functions. 
Same conventions and symbols as in Fig. [2] 



The model we explore is simplified but has the advantages 
that the results can be converged and are free of significant 
approximations. The simplified interaction used in the model 
retains essential features of the Coulomb interaction. 

We have shown (Fig. [2] and 0) that the effective densi- 
ties and potentials are explicit functions of a cost function. 
Potentials and densities very different to the exact DFT so- 
lutions are obtained if we enforce properties beyond p(r) in 
the cost function. The exact DFT wave-function matches p(r) 
with complete disregard to other elements of the many body 
wave-function structure. Since the Hohenberg-KohntheorernA 
is valid, optimizing other properties of the non-interacting 
ground state in general requires changing the potential with 
a resulting impact on the density. 

We find that mean field methods, while giving an accurate 
description of the density can mislead us in other aspects of 
the wave-function structure such as the nodal surface. In this 
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paper we argue that, among the pure-state-non-interacting v- 
representable densities, there is at least one that more accu- 
rately describes the interacting ground state nodes. We can 
optimize the wave-function associated with this density with 
a cost function. The cost function form depends on the prop- 
erty we target to retain and also the optimal density we find. 
For the nodes, the optimal cost function is clearly the DMC 
energy. A fixed cost function establishes a density-density 
correspondence which can be described as an operator U that 
transforms interacting densities into non-interacting ones. 

While finding the functional form of U dmc is a task be- 
yond the scope of this paper, we argue that we can expect this 
operator to be highly non-local and very different from the 
identity, in particular for strong electron-electron correlations. 
We find that we can improve the nodes with some simple cost 
functions, but the best nodes we found were obtained solv- 
ing the non-interacting problem without the addition of any 
effective potential. We cannot exclude the possibility that this 
result might well be an accident of the model. Our result, how- 
ever, shows that the popular expectation that the DFT solution 
is a good starting point for nodes is not valid in general. 

Optimal wave-functions can be found by altering an 
external potential. This idea is not new. In prac- 
tice wave-functions are optimized with the trial wave- 
function only minimizing the ground state energy Kvmc = 
($t| e~ J He~ J /(($t| e~ 2J |$t)). or the variance of 
the ground state energy. Replacing Kvmc into Eq. (0 
leads to a procedure similar to the optimization of Filippi and 
Fahyi^ providing additional support to that method. In the 
case of Refs. the nodes are selected by adjusting the mix 
of density functionals that gives the lowest DMC energy for 
a small system. The same mix is then used in a larger sys- 
tem. This procedure is in fact equivalent to optimizing the 
shape of the effective potential with the restriction of remain- 
ing a linear combination of two of more exchange-correlation 
potentials. We find that the change in the effective potential 
required to optimize the nodes could be of the order of the 
Hartree potential, since the wave-function with the best nodes 
is the non-interacting ground state, without any effective po- 
tential, for all the range of interaction strengths and shapes 
explored. Our results suggest that counter intuitive directions 
for potential optimizations should be explored to improve the 
nodes. 

Potential optimization has also been applied for the predic- 



tion of electronic excitations. Since p(r) determines V(r) (but 
from a constant), the excitation spectra is a function 

of p(r). This allows defining cost functions K ex to match 
the spectra of a non-interacting system. In order to minimize 
K ex one should do the derivatives 8e v / 5Vks{ v ) as m Ref 30 . 
When {E„}is taken from experiment, the search of a potential 
giving a non-interacting density that minimizes K ex is equiv- 
alent to the empirical potential method^. Unfortunately, in 
this case the electronic density can no longer be used to ob- 
tain the forces on the atoms. The existence of a single density 
functional that can be used to obtain the excitation spectra of 
any system is then a subject of debate. 

In summary, although the popular languages of electronic 
structure theory all share the same quantum mechanical un- 
derpinnings, when applied by experts to physical systems we 
often reach different conclusions. Many experts in QMC pre- 
fer HF wave-functions, while in contrast calculations done 
within the GW-BSE approach often rely on LDA derived 
wave-functions and energies 8 , while some hybrid density 
functionals obtain single particle excitations in direct agree- 
ment with excitation spectra 3 ?. We argue that as different the- 
ories need to retain different properties of the same ground 
state wave-function to minimize errors, different functionals 
should also be used. In some cases these functionals corre- 
spond to the minimization of cost functions designed to re- 
tain properties of the many body ground state in the non- 
interacting wave-function. Since some properties are favored 
at the expense of others, it is unlikely that we can use the same 
functional universally: we find that a function designed for op- 
timal nodes is a bad source of densities and vise versa. Here 
we give a qualitative picture of the size of the differences that 
one can expect as correlations start to dominate. With increas- 
ing interaction strength, the exact Kohn-Sham non-interacting 
wave-function becomes a much poorer description of several 
properties of the many-body ground state. Methods that go 
beyond DFT are limited to the nearly non-interacting limit if 
they depend strongly on DFT derived wave-functions. 
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